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O ' Abstract 

00 ' 

0^ ■ A new approach to gravitational gauge- invariant perturbation theory begins 

Q \ from the fourth-order Einstein-Ricci system, a hyperbolic formulation of gravity 

for arbitrary lapse and shift whose centerpiece is a wave equation for curvature. 
In the Minkowski and Schwarzschild backgrounds, an intertwining operator pro- 
cedure is used to separate physical gauge-invariant curvature perturbations from 
unphysical ones. In the Schwarzschild case, physical variables are found which 
^ ■ satisfy the Regge- Wheeler equation in both odd and even parity. In both cases, 

the unphysical "gauge" degrees of freedom are identified with violations of the 
linearized Hamiltonian and momentum constraints, and they are found to evolve 
among themselves as a closed subsystem. If the constraints are violated, say by 
numerical finite-differencing, this system describes the hyperbolic evolution of the 
constraint violation. It is argued that an underlying raison d'etre of causal hyper- 
bolic formulations is to make the evolution of constraint violations well-posed. 
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1 Introduction 



The paradigmatic application of perturbation theory in general relativity is to de- 
scribe distortions of the Schwarzschild black hole[jl|. Moncrief's classic gauge-invariant 
treatment [0] of this problem is a remarkable piece of work, but to many its success 
has an air of the magical to it. The principal motivation for this paper is to use the 
Schwarzschild problem to introduce a powerful and elegant new method for gravitational 
gauge-invariant perturbation theory. This method is practically algorithmic and follows 
naturally from a hyperbolic formulation of the Einstein equations^ [|, referred to 
hereafter as the Einstein-Ricci formulation, whose centerpiece is a wave equation for 
curvature. Other applications of this new hyperbolic formulation have been pursued in 
1,0- 

Besides serving as an example of a general technique, several new insights into black 
hole perturbation theory arise: The first of these is the recognition that the natural 
gauge-invariant perturbations are curvature perturbations, not metric perturbations, 
though one can get "back to the metric" if one desires. Many people have shared the 
intuition that Einstein's theory should truly be a theory of propagating curvature, but 
here that ideal is explicitly realized. Perhaps more surprising is the discovery that 
the linearized Hamiltonian and momentum constraints themselves constitute the non- 
physical gauge degrees of freedom. Violations of the constraints represent perturbations 
that take the background solution away from the constraint hypersurface defining the 
physical theory. Finally, the physical gauge-invariant variables satisfy Regge- Wheeler 
equations||, coupled to constraint violations, in both odd and even parity. 

The following schematic form is found in detail below for the system of equations 
describing the evolution of the radial part of perturbations about the Schwarzschild 
background 

ds 2 = -N 2 dt 2 + N~ 2 dr 2 + r 2 d9 2 + r 2 sin 2 6(p 2 , 

N 2 = 1 — 2M/r, on a t =const., Kij = slice. The gauge-invariant physical degrees of 
freedom are A u in odd-parity and A g in even-parity. (Their expression in terms of the ra- 
dial parts of perturbations of 3+1 variables will be given below.) Momentum constraint 
violations are encoded in c u q for odd-parity and c g g and c gr in even-parity. Violation 
of the Hamiltonian constraint is encoded in Ch- These equations are valid in a neigh- 
borhood of the Schwarzschild initial data, which is a "point" on the physical constraint 
hypersurface in the phase space of Einsteinian initial data. They reveal how pertur- 
bations of Schwarzschild behave even when the constraints are violated perturbatively. 
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where Vrw is the Regge- Wheeler potential [||. 

There are two immediate significant observations to make about this system of equa- 
tions. The constraint variables evolve among themselves under the hyperbolic sub- 
system of equations (|[), (||). If they and their time derivatives vanish initially, they are 
guaranteed to vanish for all time. This is conjectured to be a general feature of con- 
strained hyperbolic systems and is discussed further below. Next, when the constraints 
are satisfied, so the constraint variables vanish, both the odd- and even-parity physical 
perturbations satisfy the Regge- Wheeler equation ||. This somewhat surprising result 
runs counter to folklore which attributes the existence of the familiar Zerilli equation |1 
for even-parity perturbations somehow to parity. The new derivation reveals that parity 
is not the issue. 

Chandrasekhar [|l| has observed that because there is no difference between even and 
odd parity in the Newman- Penrose formalism, there is no reason to expect different 
equations for the different parities. This was part of the motivation that led him to 
construct the transformations [O] between the Zerilli and Regge- Wheeler equations and 



-d 2 + (N 2 d r ) 2 - V gt 



(2) 
(3) 



0, (4) 



their Newman- Penrose analog, the Bardeen-Press equation [|Tl|], that make their isospec- 
trality clear. Such isospectral transformations between equations are ubiquitous [O, [13 



but transformations which relate similarly simple potentials arise only under special 
conditions (cf. e.g. [jl4|), something which could well be attributed to accident. It is 



clear from the present work that the existence of the isospectral transformation allows 
one to reach both the Regge- Wheeler and Zerilli equations, in either parity, but that 
the Regge- Wheeler equation is preferable on the grounds of simplicity. In particular, 
the Regge- Wheeler equation has regular singular points only at physically significant 
locations while the Zerilli equation has a further regular singular point at an angular 
momentum dependent location. 
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An intuitive argument helps to explicate the structure of the evolution equations 
(0)"(i)- I n gauge-invariant perturbation theory, the natural candidate variables for 
perturbation are those which vanish in the background [p. 5||. In a constrained theory, 
the constraints themselves furnish a natural subset of such variables: The constraints 
are constructed from the degrees of freedom of the full unconstrained theory; they 
are satisfied, i.e. vanish, in the background; under arbitrary perturbations, they are 
generally nonzero. Indeed, just as satisfaction of the constraints defines the physical 
sector of the theory, their violation (partially) parametrizes the theory away from the 
physical sector of the theory. The variables based on the constraints will be referred to 
as constraint variables. In distinction, the variables which parametrize the constraint 
hypersurface will be referred to as physical variables. The constraint variables are a 
wise choice to describe some of the "gauge" degrees of freedom of a constrained theory, 
gauge in the sense of characterizing unphysical aspects of motion. 

By a constrained hyperbolic theory, we mean a constrained theory which is well- 
posed, and in particular one in which the constraints are guaranteed to remain satisfied 
provided they are satisfied initially. It seems intuitively clear that if such a theory were 
expressed in terms of constraint variables and physical variables, the system of equations 
would split much as do. There would be a subsystem in which the constraint 

variables evolve among themselves, and there would be further equations in which the 
physical variables couple to the constraint variables. Physical variables could appear in 
the constraint subsystem only nonlinearly, if at all, multiplied by constraint variables, 
or else their nonzero presence would act as a source to drive the constraint variables 
away from zero, given vanishing initial data. Because of this, the physical variables 
cannot appear in the constraint variable subsystem in the perturbative setting. That a 
system of equations would admit such a closed internal system of equations is obviously 
special, but it is the feature which distinguishes a constrained hyperbolic system from 
an unconstrained one. This qualitative analysis is expected to hold in the fully nonlinear 
theory, and work is in progress to demonstrate this||16||. 

One of the subtle issues that arises in constrained physical theories like general rel- 
ativity is that while the theory is hyperbolic on physical grounds, its mathematical 
representation in redundant variables may not be. In other words, were the theory re- 
duced to the true degrees of freedom, it would be a hyperbolic theory. However, when 
expressed in redundant variables, the theory may not impose hyperbolic evolution on 
unphysical combinations of variables, e.g. the constraint variables. Einstein's equa- 
tions, viewed as a system of differential equations for the metric, are an example of a 
constrained system which is physically hyperbolic but not mathematically so. If further 
restrictions are applied, e.g. through special coordinate conditions, a modified system 
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of equations can be made hyperbolic |T8 

The several new recent hyperbolic formulations follow a different route to hyperbol- 
icity|3|, [|, [§, By various extensions, the redundant sector of the theory is enlarged, so 
the whole theory becomes manifestly hyperbolic. A key virtue of achieving a hyperbolic 
formulation is that the system of equations is well-posed, so that, in particular in causal 
hyperbolic formulations, this guarantees that violation of the constraints evolves in a 
predictable, though not necessarily stable, fashion ]16||. An ill-posed theory is vulnera- 
ble to catastrophic breakdown: arbitrarily small perturbations may lead to arbitrarily 
large deviations arbitrarily quickly. In a well-posed theory, there may be exponentially 
growing modes, which may in turn be identified as instabilities, but their growth rate is 
bounded. 

In numerical simulations, the urgency of these considerations comes to the fore. 
Except in special cases, the very implementation of the constrained system will introduce 
constraint violation because finite-differencing generally does not respect the constraints. 
While there exist numerical methods for evolving ill-posed systems, there are many more 
highly developed methods for handling hyperbolic systems. One can choose to re-solve 
the constraints at each time step to try to remain near the constraint hypersurface, 
or one can face the full theory and address the issue of controlling constraint violation 
directly. In the perturbative example here, once the splitting into constraint and physical 
variables has been made, preservation of the constraints under differencing is no longer 
an issue — this reflects the fact that the theory has been reduced to the linearized true 
degrees of freedom. On the other hand, one can undertake a stability analysis of the 
system with constraint violations which will reflect on the theory in the original variables, 
and thereby gain insight into the nature of instabilities likely to appear in a numerical 
simulation of the full theory in the original variables. Work on these issues is in progress. 



2 Intertwining procedure 

We outline the procedure for gravitational gauge-invariant perturbation theory as fol- 
lows: Consider perturbations about the Schwarzschild background. The fourth-order 
form of the Einstein- Ricci formulation [H] is a hyperbolic system of a type termed "hy- 



perbolic non-strict" [ 20| . One equation of this system is a wave equation for the time- 
derivative of the extrinsic curvature, d^K^ (<% = d/dt — Cp, where Cp is the Lie deriva- 
tive along the shift (3). Because d^K^ is part of the decomposition of the Riemann tensor 
RoiOj in 3+1 variables, this wave equation propagates curvature. The time derivatives 
of the components of the extrinsic curvature vanish in the Schwarzschild background, so 
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their perturbations are necessarily gauge- invariant. They are chosen as the perturbative 
quantities and will be referred to hereafter as curvature perturbations. Their selection 
accords with the principle that the natural candidates for gauge-invariant quantities are 



those which vanish in the background ||15||. Note, in contrast, that the Riemann tensor 
itself -RoiOj does not vanish in the background. The Einstein-Ricci system is hyperbolic 
for arbitrarily specified lapse and shift, which consequently do not have to be perturbed. 

There are six components of doKij, and one has six coupled wave equations in terms 
of them. In addition, there are four "constraints," formed from time-derivatives of the 
three momentum constraints and the Hamiltonian constraint. This leaves two indepen- 
dent equations to be found. If one makes a tensor spherical harmonic decomposition of 
the perturbations of doKij, the latter naturally divide into odd and even parity. The 
two dynamical equations and one constraint with odd parity lead to the Regge- Wheeler 
equation. This leaves four dynamical equations and three constraints with even parity. 

The isolation of a wave equation from the even parity system is the technical crux of 
the gravitational perturbation calculation. An analogy will make the procedure clear. 
Suppose that one had four linear algebraic equations in four unknowns. To solve them, 
one would take linear combinations of the equations to isolate linear combinations of 
the unknowns. Equivalently, one diagonalizes the matrix of coefficients in the equations 
by transforming to a different basis. 

In the present case, we have four linear differential equations in four unknowns. We 
construct linear differential combinations of the equations to isolate linear differential 
combinations of the unknowns. We do not diagonalize the matrix of differential op- 
erators, but we re-group the variables to split the differential equations into partially 
uncoupled form which is diagonal when the constraints hold. We do this by grouping 
the original variables into combinations which constitute the constraints and new vari- 
ables which are independent of the constraints. This is accomplished by an intertwining 
transformation [OL [13| in which one matrix of differential operators M\ is transformed 



into another M 2 by a matrix operator D 

M 2 V = VM X . (5) 

If T> were invertible, one would have the familiar expression T>M{D~ l = M 2 for a basis 
change. When the constraints hold, the matrix M 2 is diagonal; when they do not, the 
matrix M 2 is simply a matrix in a different basis. 
Nonperturbatively, the momentum constraints are 

Roj = -N(w k K ]k - VjH) = (6) 
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where H = K \ is the trace of the extrinsic curvature and overbars indicate spatial 
quantities, here spatial covariant derivatives. The perturbed time-derivative of the mo- 
mentum constraints are differential linear combinations of the curvature perturbations 
doKij. They are thus appropriate combinations to use in re-organizing the system of 
four linear differential equations. As displayed in (|3|), (f|), the momentum constraint 
variables evolve among themselves together with the Hamiltonian constraint variable. 
The Hamiltonian constraint itself 

G% = -^(H 2 -K i ^K lJ +R) = 

involves only metric perturbations when perturbed about a = slice (where the 
underline indicates a background quantity). As a spatial constraint, this constraint may 
be violated when evaluated on a general metric and extrinsic curvature. This reflects 
violation of the Hamiltonian constraint. Taking two time derivatives converts the metric 
perturbations to curvature perturbations and gives an equation for the evolution of 
violation of the Hamiltonian constraint. This is simply the time derivative of the doubly- 
contracted Bianchi identity V^G^o — 0, which can therefore be read as an equation 
evolving violations of the Hamiltonian constraint through coupling to violations of the 
momentum constraints (cf. pHf). 
Separately, the relation 

R° = -N~ l d Q H + KijK ij - N- l V k V k N = 

stands as a dynamical "constraint" coupling curvature and metric perturbations in per- 
turbation theory. This equation is not part of the Einstein-Ricci system and hence may 
be violated. Two time derivatives again convert the metric perturbations to curvature 
perturbations. 
The identity 

R k k = — 2G°o + -R°o (7) 

relates the two constraints to the field equations. The second time derivative of R k k 
occurs as part of the dynamical equations which define the fourth-order Einstein-Ricci 
theory and is therefore part of an equation which does not admit violations. Thus, two 
time derivatives of (0) establishes the correlation between violations of G°o = and of 
R°o = 0. In view of this correlation, through a useful abuse of language, R°q will be 
referred to also as the Hamiltonian constraint, and the so-called Hamiltonian constraint 
variable in (|l])-(fD reflects violations of the R°q constraint. Violations of £7% can be 
eliminated from the second time derivative of (^) using the doubly-contracted Bianchi 
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identity mentioned above, resulting in the wave equation in (ffl) for violations of R°q = 0. 
Details of this will be given below. 

Thus, three of the four equations have been "diagonalized" already, and it is not 
necessary to begin from a general matrix ansatz for the intertwining matrix T>, though 
this could have been done. Let g k (r, t), k — 1, . . . , 4, be the radial parts of the even-parity 
curvature perturbations, and let L k (g 1 , . . . , g±, r, t) = 0, k = 1, . . . , 4, be the associated 
dynamical coupled second-order differential equations for them. (Note that L k is not 
simply the differential operator, but the full equation; this is indicated notationally by 
the hat.) It is sufficient to consider a first-order differential linear combination of the four 
dynamical equations to achieve a single wave operator acting on the same differential 
linear combination of radial curvature perturbations 

4 4 

^(b k d r + a k )L k = L RW ^{b k d r + a k )g k + /(-R o, Roj), (8) 

k=l k=l 



where 

Lrw 



-d 2 t + (1 )d r {\ )d r - V RW {r) (9) 



2M 2M, 
— )d r (l 

is a wave operator of Regge- Wheeler form and / is a linear differential operator applied 
to the constraints. The potential 

V RW r = 1 - — 10 

is the Regge- Wheeler potential — it could be left undetermined initially and solved for 
self-consistently. (Alternatively, one could attempt to reach the Zerilli potential starting 
from a second-order differential combination.) The a, and b{ are functions of r to be 
determined by equating like coefficients of derivatives of the g k s term-by-term. 

A simplification is achieved by recognizing that the momentum constraints can be 
used to eliminate linear spatial derivatives of particular radial coefficients in favor of 
other terms. This allows one to set two of the b k to zero without loss of generality. The 
equations for the a k and b k obtained from (^) by equating terms are over-determined. 
After determining a k and b k from a sufficient set of equations, the remaining equations 
are consistency conditions. If the equations are not consistent, an intertwining of the 
chosen form does not exist, though one of a different form, e.g. higher order differential 
or infinite order by means of integral operators, may exist. 

The operator transformation between the Regge- Wheeler and Zerilli equations that 



Chandrasekhar|l|] discusses is also an example of intertwining fll4|, |22]| . Because inter- 



twining transformations compose, when diagonalizing the system above, either equation 
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can be reached by using an appropriate linear differential combination of the curvature 
perturbations. One has 

D 1 L RW = L Z D 1 , (11) 
D 2 L Z = L RW D 2 , 

where the Regge- Wheeler operator L^w is given by @ and the Zerilli operator is 

L z = -d 2 + (1 - ™)d r (l - ^-)d r - V z {r) (12) 

with the Zerilli potential [1], |J 

2N 2 ( \ 

V z (r) = — (n + l)n 2 r 3 + 3MnV + 9M 2 nr + 9M 3 , (13) 

r 6 {nr + 3M) 2 \ J 

where N 2 = 1 — 2M/r and n = (£ — 1)(£ + 2)/2. It is important to note that one can 
pass in either direction by means of differential operators; the inverse transformation 



need not be an integral operator. The intertwining operators are [14 



2ilf 3M(r-2M) . , 

Di = 1 )d r + K ( +u>, 14 

r r z (nr + SM) 

2M 3M(r - 2M) 

° 2 = 1 ^ ~ 2/ —^TK 

r r z [nr + SM ) 

where uo = n{n + l)/(3Af) = - 1)£{£ + 1){£ + 2)/(12M). The intertwining relations 
(O) become more transparent when one recognizes that they are a consequence of 
associativity, Di(D 2 Di) = (DiD 2 )Di and D 2 (D 1 D 2 ) = (D 2 D 1 )D 2 , because the radial 
parts of the Regge- Wheeler and Zerilli operators factorize 

L RW = D^-d^+u 2 , (15) 
L z = D 1 D 2 -d 2 +tu 2 . 

This factorization property is very suggestive, but it turns out not to be the best property 
to generalize intertwining. 

3 Einstein-Ricci formulation 

As shown in P, f|, [|] , the dynamical part of Einstein's equations can be cast in hyperbolic 
form in 3+1 language as the Einstein-Ricci formulation. Consider a globally hyperbolic 
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manifold of topology £ x R with the metric 

ds 2 = -N 2 dt 2 + g ij (dx i + ftdt) (dx j + f3 j dt), (16) 

where N is the lapse, p 1 is the shift, and is the spatial metric. Introduce the non- 
coordinate co-frame, 

e° = dt, 6 l = dx l + ftdt. (17) 
with corresponding dual (convective) derivatives 

d Q = d/dt - pd/dx\ di = d/dx\ (18) 

The natural time derivative for evolution is [EH 

4 = d + (3 k d k -c p = d/dt - C p , (19) 

where Cp is the Lie derivative in a time slice S along the shift vector. In combination 
with the lapse as N~ 1 8q, this is the derivative with respect to proper time along the 
normal to S. The extrinsic curvature Ky of E is defined as 

B o9ij = -2NKij. (20) 

From this follows the equation for the evolution of the Christoffel connection 

dof \ 3 = -g mk [V 3 {NK im ) + Vi(NK mj ) - V m {NK l3 )]. (21) 

Barred quantities are three-dimensional. 

The dynamical Einstein equations, Rij = pij, where pij is a matter source, are 
equivalent to a third-order Einstein- Ricci system. The system is labelled "third-order" 
by the equivalent of the highest number of derivatives of gij that can occur in the theory. 
(These derivatives need not and generally do not appear explicitly.) The third-order 
Einstein- Ricci system is obtained from (20) and a wave equation for K~ij 

N 2 uK i3 + NJ i3 + NSij = NSlij = NidoRij - ViR oj - V 3 R 0i ), (22) 

where 

□ = -N^doN^Bo + V fc V fc . (23) 

Here, Qij is a matter source term which, as a consequence of the field equations, vanishes 
in vacuum. For this paper, we restrict attention to vacuum spacetimes. is a nonlinear 
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self- interaction term. If we denote the trace of the extrinsic curvature by H = K \, Jij 
is given by 

= d (HK ij -2K i k K jk ) + (N-%N + H)V i V j N 
-2iV- 1 (V fe iV)V( i (iV^ i) ) + 3{V k N)V k K l:j (24) 
+N- 1 K ij V k (NV k N) - 2V {i (K j) k V k N) + N^HV.VjN 2 
+2iV- 1 (V (j #)(V ?) iV 2 ) - 2NK k {l R j)k - 2NR kijm K k 



where M(y) = + Mji). The three-curvatures can be expressed in terms of four- 

curvatures and then eliminated using the field equations. 
Finally, Sij is a slicing term, 

S {j = -N^ViVjidoN + N 2 H). (25) 

This must be equal to a functional involving fewer than third derivatives of the metric 
to assure the hyperbolic (wave) nature of the equation (p2|). Two simple ways to do this 
are the following. One may impose the harmonic slicing condition 

d N + N 2 H = 0, (26) 

in which case 

= 0. (27) 

Alternatively, one may specify the mean curvature by demanding H = h(x,t). The 
lapse function iV is then determined by solving the elliptic equation 

V fc V fc iV = -d Q h(x, t) + NKijK ij . (28) 

The special case of maximal slicing, H = 0, gives 

S { j = -N'^iVjdoN. (29) 

In this case, the Einstein- Ricci system is mixed hyperbolic-elliptic [p4 



With harmonic slicing, the dynamical part of Einstein's equations are given by the 
definition of the extrinsic curvature (p0|) , the wave equation (22) , and the harmonic 
slicing condition (|26|). (In this case, all the equations of motion are equivalent to a 
first-order symmetric hyperbolic system[|3], |J with characteristics the light cone and the 
direction orthogonal to the time slices.) With slicing given by specified mean curvature, 
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the dynamical equations are fl2"0|) and (0), and the lapse is determined by interleaving 
the solution of the elliptic equation (^) on each time-slice. 

To complete the new formulation of the Einstein equations, we specify the initial 
Cauchy data. Initial and K v - must be chosen compatible with the Gauss-Codazzi 
(Hamiltonian and momentum) constraints. The Hamiltonian constraint is 

G°o = l(R°o - R\) = -\{R + H 2 - K mk K mk ) = 0, (30) 
and the momentum constraints are 

NG° l = V k (K k l -5 k i H) = 0. (31) 



These are treated as an elliptic system on the initial slice by the usual methods |23"1, |25 
Furthermore, to guarantee that the dynamical equations produce a solution to Einstein's 
equations, it is also necessary that d§Ky be specified initially so that Einstein's equations 
hold on the initial slice 

R i3 = R i0 - N^doKij + HKij - 2K tk K k J - N' 1 V.V.iV, (32) 

where Ry on the left hand side is replaced by its expression in terms of the matter 
sources. In the case of harmonic slicing, the lapse is specified on the initial slice. The 
shift, hidden in do, is a freely specifiable function on spacetime and is not a dynamical 
variable of the system. 

To obviate the need for special handling of the slicing term and to allow the lapse 
to be specified freely, one can take another time-derivative, apply a further constraint, 
and pass to the fourth-order Einstein-Ricci formulation. These equations proceed from 
(p0|) and the equation 

Ndo(NnKij) + NdoiJij + Sij) + NViVjiNdoH - N 2 K mk K mk + NV k V k N) = 

= N&ij = N0lR l3 - 2d V (i R j)0 + ViVjRoo). (33) 

The effect of including i?oo is to incorporate the Hamiltonian constraint while cancelling 
the threatening highest derivative term of H in S^. The Cauchy data are extended by 
specifying d^K^, subject to the requirement (in vacuum) that doRy = hold on the 
initial slice. Both the lapse and the shift are freely specifiable functions on spacetime 
and are not dynamical variables. The system (p0|), ( |33"D is hyperbolic non-strict in the 
sense of Leray-Ohya|[20||. In particular, it is well-posed with solutions in an appropriate 
Gevrey (not Sobolev) class[H]. 
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4 Linearized gravity in the hyperbolic formulation 



Let us first consider a weak-field analysis of the equations above around a flat spacetime 
background with a Minkowski metric. Let gij = g.^ + g'^ and Kij = + K'ij where 
underlines denote background values and primes the first order corrections. In the 
present case, g = 5ij and = 0. In the fourth-order Einstein-Ricci formulation, the 
lapse and shift can be arbitrarily specified: allow a lapse perturbation N = N + N' = 
1+N', but set the shift perturbation to zero, f3 % = (3 l = 0. The use of a lapse perturbation 
is unnecessary, but its inclusion will demonstrate this explicitly. (The lapse perturbation 
N 1 here should not be confused with N' = d r N used later in the paper.) At first order, 
( pOD becomes 

d t g' ij = -2K' ir (34) 



The wave equation (|33|) from the fourth-order Einstein-Ricci formulation is 

ndtK'ij - Vi % (d 2 t N' + d t H') + % VjR'oo = (35) 

(o = —d 2 + V. fc V fc ), where V, is the spatial covariant derivative with respect to the 
background metric q ... 

The weak-field form of the identity 

Rm = N 9oKij + N 2 K ik K k , + NViVjN (36) 

is 

& m = d t K' ij + V& j N'. (37) 

Its trace gives 

R'w = d t H' + V_ k V k N'. (38) 

Using this in (|35D gives 



n(d t K'ij + YiVjN') = m' m = 0. (39) 

This equation is clearly lapse independent. Furthermore, the perturbations are explicitly 
Riemann tensor perturbations. 

There are six degrees of freedom to R'cuoj, but we expect only two physical gauge- 
invariant degrees of freedom. The resolution of this apparent paradox is to recognize 
that by tracing ( |39"|) one obtains an equation for R' 00 

□#oo = 0, (40) 
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while taking a (spacetime) divergence and using the contracted Bianchi identity, 

V J -R iOj = — Vo-Roj + Vj-Roo, (41) 
one reaches an equation for d t R'oj 

nd t R' 0j = 0. (42) 

Thus, four of the six equations evolve constraints. A general perturbation may violate 
the constraints, and these four equations determine the constraint violation evolution. 
The physical degrees of freedom are the remaining two degrees of freedom of R'oioj, the 
transverse traceless parts. 

To make contact with traditional analyses, one can examine (|37| ) in more detail by 
reading it as an equation for dtK'ij and splitting K'ij into a sum of transverse traceless, 
longitudinal traceless and trace parts 

K% = K'JJ + K'f + l - 9ij H>. (43) 

The trace of fl37D is 

d t H' = R' 00 -V k V k N' (44) 

as above. Note that R'oo is not set to zero as traditionally done. This is because a 
general perturbation will violate the constraint, and it is easier to work with free rather 
than constrained perturbations. 

The longitudinal part can be identified by using the definition that the divergence of 
the transverse part vanishes. Split the extrinsic curvature perturbation and the Riemann 
tensor into a sum of transverse traceless and longitudinal (with trace) parts 

K'^ = K'JJ + K' i: j, R'oioj = R'oioj + R'oioj i ( 45 ) 

where the divergence of the TT parts is assumed to vanish. The divergence of (|3~7| ) is 
then 

d&K% = V { R% oj - YX V, .V. (46) 
Using the definition of the perturbative momentum constraint 

&oj=-&K' ij -%H') (47) 

and ( JH| ) leads to the perturbative form of the contracted Bianchi identity (P"T|) 

— d t R'oj = V^-R'oioj — Vj-R'oo- (48) 
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This reveals the longitudinal nature of this identity. Stripping the divergence away from 
(|46|) and removing the trace gives the tracefree longitudinal equation 

d t Kf = R' h mj - i^'oo - (V^iV' - l -g l3 V k V k N'). (49) 

Finally, the transverse traceless equation is 

d t K'1J = R%. (50) 
If one further splits ([34]) to obtain the transverse traceless part, then one finds 

- \d?gf = R'Zj, (51) 

which is a well-known result. One of the virtues of this analysis is that it clarifies the 
role of the transverse and longitudinal parts of the Riemann tensor and emphasizes that 
the true physical degrees of the freedom of the linearized gravitational field lie in the 
transverse traceless part. 



5 Gauge- invariant perturbation theory: Schwarzschild 

The fourth-order Einstein-Ricci system is the natural one for gauge-invariant perturba- 
tion theory. Because the lapse and the shift are freely specifiable, their perturbations 
do not need to be considered. Furthermore, because one is working with perturbations 
which have the dimensions of curvature, it is easier to find variables which are gauge- 
invariant in the background, and hence natural candidates to perturb. It is convenient 
nevertheless to begin by perturbing the third-order Einstein-Ricci system (without fixing 
the slicing term) and later to take a time derivative to reach the perturbed fourth-order 
theory. This procedure organizes and simplifies the computation because the background 
is time-independent. 

As perturbations, let = g.. + g'^, Kij = + K'ij, T l jk = jk + f n jk, where an 
underline indicates the background quantity and a prime the first order perturbation. 
In the fourth-order theory, the lapse is unperturbed, so iV = N_, and the underline will 
be suppressed. The shift (3 k = (3 k is also unperturbed. 

To be explicit, consider perturbations of the Schwarzschild background 

ds * = -(1 - —)dt 2 + (1 - — yV + r 2 d6 2 + r 2 sin 2 #0 2 . (52) 
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This example admits three useful simplifications: The background shift vanishes, (3 k = 0, 
so do = dt- The background lapse N is time-independent, so many terms in ( p^ ) and 
(|53]) vanish. Finally, the background extrinsic curvature and its derivatives (space 
and time) vanish on the natural slice (t =constant) to be perturbed. This implies that 
their Lie derivatives along an arbitrary vector v vanish, e.g. 

CvKn = v k d k K l3 + K kj d t v k + K ik d jV k = 0. (53) 

Thus, they are perturbatively gauge-invariant variables, and therefore natural candidates 
for perturbation. Their selection accords with the general principle emphasized in ]15| 



that one should always choose to perturb quantities which vanish in the background. 
In addition, many more terms in and ( |3"B| ) vanish because the extrinsic curvature 
is zero in the background. 

The result of perturbing fl22|) in a background sharing these properties (time-inde- 
pendent lapse, vanishing shift, and vanishing background extrinsic curvature) with 
Schwarzschild is the wave equation 

M^. = -dlK'v + N 2 V k V k K' l3 - N 2 V_^H' - ANV_ k N% t K' j)k (54) 
+3N(V k N)V k K' ij - 2NK' k{i V j) V k N - 2K' k{i (% ) N) V_ k N 
+NH%V j N + K'ij(V_ k N)V_ k N - 2N 2 K' k(i R j) k - 2N 2 R k ij m K' km . 

(An additional term K' i jNv!'V_ k N vanishes in the Schwarzschild background and is 
not included.) The spatial Ricci and Riemann tensors are computed in the background 
metric. 

Note the presence of the slicing term as the second derivative of H' = trK'. This 
spoils the hyperbolicity of the perturbative wave equation. If one wished, one could 
return to the full theory and impose harmonic slicing. When perturbing the theory, 
one would have to allow perturbations of the lapse, but this second derivative term 
would be removed. Interestingly, the perturbed lapse would still not appear in the 
perturbed wave equation. If one went on to introduce a tensor spherical harmonic 
multipole decomposition of the perturbations of the extrinsic curvature, as will be done 
below for the time-derivative of the extrinsic curvature, perturbed third-order equations 
would result. These form the basis of the perturbative- matching outer-boundary module 
being used as one method of extracting gravitational waves and imposing outer boundary 
conditions on the three-dimensional numerical simulations carried out by the Binary 



Black Hole Alliance [26 . 
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The second derivative of H in the full theory can also be removed by passing to the 
fourth-order theory by taking a time-derivative and adding second spatial derivatives of 
R 00 . Because the background is time-independent and the background extrinsic curva- 
ture vanishes, the effect of taking a time-derivative is simply to replace the perturbed 
extrinsic curvature K'ij by its time derivative d t K' \j = K[-. The additional perturbed 
i?oo term has the form 

NVS-jR'oo = N 2 V_S-jH' + NH'VjVjN + 2N(V (l N) V j} H' (55) 

-NVS. j (Ng mn g' nj ^ k ZS*N + Ng^f'^diN) 

One sees that by using -R'oo, the second derivative of H' is cancelled, but as well 
metric perturbations are introduced into the perturbed fourth-order equations through 
g' and f '. This is the price that one has paid to achieve explicit hyperbolicity in the 
full theory. The next step is to eliminate the metric perturbations in favor of time- 
differentiated extrinsic curvature perturbations. This is done by using -R'oo- On first 
sight, this appears simply to undo the step which cancelled the second derivative of H, 
but this is not so. We desire to understand the behavior of the theory away from the 
constraint hypersurface, so we do not set -R'oo = 0, but leave it as a free variable which 
happens to vanish on-shell. The presence of this off-shell term marks the fact that the 
second derivative of H has been cancelled in the full theory. 

This raises a valuable point. It is important to emphasize that one starts from a 
well-posed hyperbolic theory before perturbing. If the Rq term were not present in 
the full theory, as it is not when one simply takes a time-derivative of the third order 
theory (without fixing the slicing), one would have a system which is believed not to 
be hyperbolic. The perturbed form of that non-hyperbolic system would naively agree 
with the one we have just found if one were to set -R'oo = 0. This perturbative theory 
therefore agrees with general relativity on-shell, but the Hamiltonian constraint has not 
been fully incorporated, so that when violations of the Hamiltonian constraint occur, 
the theory cannot respond properly, that is, in a manifestly hyperbolic fashion. 

The i?oo term compensates behavior of the second spatial derivative of H to make 
the theory hyperbolic. When i? o vanishes, this compensation is evidently insignificant. 
As the magnitude of the second spatial derivatives of -Roo increase however, one suspects 
that the compensation becomes more important. Because this condition involves second 
derivatives, small variations in R 00 from zero can nevertheless produce arbitrarily large 
contributions to the equation. This reveals an important caveat: it is dangerous to 
study a constrained theory solely within the constraint hypersurface and not to consider 
its behavior when the constraints are violated. The full mathematical character of a 
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fundamental physical theory with constraints, in particular its hyperbolicity and well- 
posedness, involves off-shell information. 

Having given this warning, in the case of perturbations of Schwarzschild at least, the 
role of the R 00 term appears to be relatively innocuous. The system of equations does 
not change character dramatically when this term is removed. This suggests there is 
further room to explore the role of the i?oo term. 

The equation for the perturbed fourth-order theory is 

N% = -dlK'^ + N 2 V fc V fe ^, - iV 2 V. 4 V^' + NVS-jR'oo (56) 
-4:NV k NV {l K' j)k + 3iVV*W fc ^ - 2Nk' k[j V j) V k N + NH'V^N 
-2^ (j (V, ) iV)V /c iV + K'^N^N - 2N 2 K' k{i R j) k - 2N 2 R\ 3 m k' km . 

(Again, an additional term K'^NV V^iV vanishes in Schwarzschild case and has been 
dropped.) 



6 Odd-parity perturbations 

Because the wave operator on the Schwarzschild background separates in spherical coor- 
dinates, perturbations in the Schwarzshild background are naturally handled by making 
a multipole decomposition in tensor spherical harmonics. These in turn naturally di- 
vide into odd and even parity. Following Moncrief , we decompose the odd-parity 
perturbations as 

k'y = ux{t,r){e x )ij + u 2 {t,r)(e 2 )ij, (57) 
where Ui(t,r) and U2(t,r) are radial perturbations and 
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are the odd-parity tensor spherical harmonics. In these formulae, Yg m (9,(/)) are the 
standard scalar spherical harmonics satisfying 

(d 2 + cot Ode + — = -1(1 + l)Ytm. 
sin 9 

The notation "symm" indicates the matrices are symmetric. 

From fl56D one obtains equations for u\ and u 2 . To simplify the expressions it is 
useful to introduce 

\i = £(£ + l) 

as the (negative of the) eigenvalue of the spherical harmonics and to replace second and 
higher derivatives of N by equivalent expressions in terms of N and N' = d r N (not to 
confuse the prime with the perturbative part). From the r#-component of fl56|), after 
dividing out the common angular factor —d^Y^ m / sin 9, one has 

N 2 



d 2 tUl + iV 4 <9 2 «i + AN 3 N'd rUl + —(-Xt - 4iV 2 + 6rJVJV' + r 2 N' 2 ) Ul (60) 



(2 - X e )N 2 



+ - u 2 = 0. 

From the 6(p equation, after dividing out 

(- sin 9d%Y im + cos 9d Y im + esc 9d^Y im )/2, 

one has 

27V 3 

- d 2 t u 2 + N 4 d 2 u 2 + {-N + 2rN')d r u 2 (61) 

TV 2 47V 3 
+ — ( y -\ e + 4N 2 -4rNN' + r 2 N' 2 )u 2 + —( y -N + rN')u 1 = 0. 

Of the constraints, only one is present in odd parity. It can be found from R'oe, and 
its radial expression is 

sin. 9 R'oe „„ ,2N 3 9 , N \g - 2 

— — ™ = c ue = N 3 d rUl + {— + N 2 N') Ul + ^Nu 2 . (62) 

(A degenerate expression is found from R'q^.) If c u g = 0, the odd-parity part of the 
momentum constraints hold; if not, c u g measures the violation. The variable u 2 can be 
eliminated from ( pD| ) using this constraint, and one has 

2N 4 N 2 2N 
-d 2 u 1 +N i d 2 u 1 + UN 3 N'+ )d r u 1 + — (-\ £ +8rNN'+r 2 N' 2 )u 1 c u9 = 0. (63) 
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By rescaling u\ } the differential part of the operator can be brought into the familiar 
Regge- Wheeler form, 

-d 2 + N 2 d r N 2 d r . 

The new variable 

A u = rN Ul (64) 
is the gauge-invariant variable which satisfies the Regge- Wheeler equation 

N 2 

The familiar Regge- Wheeler potential is explicitly given by 



-d 2 + N 2 d r N 2 d r -—(Xe- QrNN') )A U - 2N 2 c uB = 0. (65) 



y Rw(r) = ^ _ m = (1 _ ^)(fcl> _ (66) 

The c u e term reflects how the Regge- Wheeler equation is modified when the odd-parity 
momentum constraint is violated. 

From the perturbative form of (|36| ) in the Schwarzschild background, the gauge- 
invariant variable A u is readily identified as the radial part of one of the components of 
the Riemann tensor 

sin OR'oroe 

A u = rNui = -r (67) 

The "extra" factor of r is an artifact which arises because the differential part of the 
Regge- Wheeler operator is written as if space were one-dimensional whereas the operator 
actually comes from the radial part of a three-dimensional Laplacian. The factor of r is 
just the scaling to change the apparent dimension of the background space. 

The equation which evolves c u q is found by forming the differential linear combination 
of ( |6"0"D and (|6T| ) implied by ( |62"D to create a second time derivative of c u q. One finds 

/V 2 A \ 

-d 2 + N 2 d r N 2 d r - — p J c u9 = 0. (68) 

(The scaling of c u g was chosen earlier to put the differential part of the operator in 
Regge- Wheeler form.) 

The odd-parity violation of the momentum constraints propagates hyperbolically on 
its own, and the pair (|65| ) and ( |S8|) together form a hyperbolic system. If c u e = and 
dtC u o = initially, ( pq) guarantees that they remain so. It is a necessary feature that the 
constraints evolve amongst themselves without involving the physical degree of freedom. 
This enables zero initial data to remain zero. We will see this again in the even-parity 
case. 
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7 Even-parity perturbations 

Following Moncrief, but relabelling the radial variables and rescaling the g\ variable to 
have more rational dimensions, we can decompose the even-parity perturbations as 



K = rg 1 (t,r)(f 1 ) ij + g 2 {t 1 r)N- 2 {f 2 ) ij + r 2 g 3 (t,r)(f 3 ) i:j + r 2 g,{f A ) %v (69) 



where 



/ dgYi m d^Yim \ 
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(73) 



\0 an 2 9Y em J 
( 

dgY tm (dgd^ - cot 9d^)Y tm 
\0 symm (9| + sin 8 cos 9de)Yi m J 
are the even-parity tensor spherical harmonics. The trace of is found to be 

H' = {g 2 + 2^3 - >^m)Yem- 
Finally, let 

The fourth-order Einstein-Ricci equations for even-parity radial perturbations 
have the following forms. From the r#-equation, up to a factor of rdgY£ m , the equation 
for gi is 



(74) 



L x : -£ gi + N*d 2 r9l + 



,2N 



N 2 



+ AN N')d r gi - —d r (g 2 + 2g 3 - \ m ) 
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N 2 

+—( y -X e -4N 2 + 10rNN' + r 2 N' 2 )g 1 (75) 

iV A^ 2 iV AT 

+ -(3iV - 2rN')g 2 + — (A, - 2)g A + -d r c h - -c h = 0. 

From the rr-equation, up to a factor of Ye m /N 2 , the equation for g 2 is 

27V 4 

L 2 . _ jV 4 ^(2(/3 - Atf 4 ) + ( iV 3 iV')«9^ 2 - iV 3 iV'^(2^3 - A^ 4 ) 

r 

+^T 1 9i + -^(-^ ~ 4iV 2 + QrNN' - r 2 N' 2 )g 2 (76) 
2N 3 



+^(N - 2rN')(2g 3 - \ m ) + iV 3 ^ + N 2 N'd r c h = 0. 
From the 00-equation, up to a factor of r 2 (ded^Yem — cot 9d^Yg m ), the equation for g 4 is 



2N 3 T „ „„„ T „, x N 2 2AT 



2 



L 4 : -9^ 4 + iV 4 «9 2 s 4 + (iV + 2rN')d r g, + —(N - rN') 9l - —g 2 - —g 3 

ry /y^i ry£ ry£ 

N 2 N 
+ — {2N 2 + ArNN' + r 2 AT ,2 )^ 4 + —c h = 0. (77) 

ry A ry £k 

Finally, taking the trace of the 99- and 00-equations, dividing out the common factor 
2Y£ m , and adding A^/2 times ([H]), one obtains the equation for g 3 

/V 4 A 7V 4 

L 3 : -d 2 g 3 + N 4 d 2 g 3 - —d r g 2 + AN*N'd r g 3 + ^-d r g A (78) 

/V 3 N 2 
+ — (2A^ - 5rN')g 2 + —(—\ t - 2N 2 + QrNN' + r 2 N' 2 )g 3 

ry Zi ryA 

A /V 3 /V 3 
+ ^ r {2N - rN')g A + — d r c h = 0. 

The (time-differentiated) perturbed momentum constraints in even-parity reduce to 
the two expressions 

R' 1 

"' -c r = d r (2g 3 -X e g A ) + -(Xtg 1 -2g 2 + 2g 3 -X e g 4 ), (79) 



NY lm ' v ^ r 

n ' w = c g6 = -rN 2 d r9l - 9l (3N 2 + rNN') + g 2 + g 3 - g 4 . (80) 



Nd e Y im 
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(Again, a degenerate expression is found for R' ^.) These define variables c r and c g g which 
are useful for studying the effects of momentum constraint violations in perturbation 
theory. By forming the appropriate differential linear combinations of (|75|)-(|78D, one 
finds that the momentum constraint variables satisfy the following system of coupled 
wave equations 



L 2N 3 ,„ . _ AN 2 N' 



- d 2 c r + N 4 d 2 c r + (N + ArN')d r c r + d r c h (81) 



r 

N 2 , „ , 9N 2A^iV 2 X f N' 

+ — (-2-\ e + 9r 2 N' 2 )d r + ^^c 9 e-^c h = 0, 



<9 2 c, e + N 4 d 2 c g9 + W 3 N'd r c ge - N 2 N'd r c h (82) 

N 2 27V 3 
+ — (-X e -2rNN' + r 2 N' 2 )d g e + — -(N + rN')d r - ,! 

Introducing rescaled momentum constraint variables 



r N 2 R' 

Cr = rN 3 c r = (83) 

and 

c 9 e = N~c ge = -^L_ (84) 
puts the differential operator part in Regge- Wheeler form, leaving the equations 

N 2 \ 
-d 2 + N 2 d r N 2 d r - — (2 + A^ + 2rNN') J c r (85) 



and 



+ ^Cge + 4N 5 N'd r c h - c h 



d 2 + N 2 d r N 2 d r - -^-)c g e + ^(N + rN')c r - N 3 N'd r c h = 0. (86) 



An evolution equation for the Rq Q constraint is found by taking the trace of 
The perturbed constraint is 



Y em = R'oo = NH' - Ng mn g ' n • g j fc V m V fc N - Ng mk t H mk %N. (87) 
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This involves metric perturbations through the g' and T' terms. A time derivative of T' 
only involves perturbed extrinsic curvatures and a second time derivative involves K[p 

d 2 r n mk = - £ « (v m (Nk; k ) + v k {Nk' ]m ) - %{Nk' mk )) . (88) 

Taking two time derivatives of fl87D allows second time derivatives of H' to be replaced 
by second derivatives of Ch and leads to the equation 

+ + + NV)3 r ~^)c h+ ^ = 0. (89) 

An alternative derivation which reveals the role of the Hamiltonian constraint G°o 
is given as follows. Let 

ChaYtm = G'\ = (90) 

Two time derivatives gives 

dUc hG Y, m ) = -V k R' k0 - N-\V k N)R' k0 , (91) 

which is simply the time derivative of the contracted Bianchi identity V^G^o as men- 
tioned above. The trace of ([56]) corresponds through ( p3|) to 

d 2 t R'\ - 2V_ k R' m + V h V k R 00 = 0. (92) 
The second time derivative of the perturbed form of (0) gives 

d 2 t R'\ = -N- 2 d 2 (c h Y em ) - 2d 2 (c hG Y em ). (93) 
Substituting this and (|90|) into (^) gives 

- N- 2 d 2 (c h Y em ) + V_ k V k {c h Y lm ) + 2iV- 1 (V fc iV) J R' fc o = 0, (94) 

which is readily confirmed to agree with (|S9D. 

The constraint variables c r , c g o, and together form a hyperbolic system. If they 
and their first time derivatives all vanish, the initial data for this system is identically 
zero and it remains zero under evolution. If any are non-zero, the constraint violations 
evolve hyperbolically. Below we will see that c r , c g o and Ch occur as source terms in the 
even-parity Regge- Wheeler equation. 



24 



8 Intertwining and the even-parity Regge- Wheeler 
equation 



Having found how constraint violations evolve, it remains to find how the physical even- 
parity degree of freedom propagates. The procedure is to form an arbitrary first-order 
differential linear combination of ([T5|)-([78|) and require that it equal a Regge- Wheeler 
wave operator 

a RW = -d 2 + N 2 d r N 2 d r -V RW (r) (95) 

acting on differential linear combination of the g k , plus constraint variables. More 
precisely, one requires 

Y,(h{r)d r + a k (r) )L k = U RW ( ^ b k (r)d r g k + a k (r)g k ) + f(c r ,c gd ,c h ). (96) 

As discussed in the introduction, this equation can be understood as a component of 
a matrix intertwining relation between a matrix differential operator in the basis of 
the g k and one in the basis which isolates the physical degree of freedom A g from the 
constraint variables c r , c g e and Ch- The other components of the matrix transformation 
were determined above when the appropriate combinations of the g k were found which 
produce equations involving only the constraint variables. 

One can begin with the potential V R w{r) of the wave operator undetermined and 
solve for it self-consistently, but this complicates some of the intermediate expressions, 
so the ansatz is made that the potential has the Regge- Wheeler form (|B6l) . The ansatz 
for the transformation (|96|) can be simplified slightly before beginning. Because the 
momentum constraint variables are built from radial derivatives of g\ and 2g 3 — \gg^ 
one can without loss of generality set &i(r) = and 63 (r) = because the momentum 
constraint variables can be used to eliminate the radial derivatives of g\ and g%. 

The procedure is now to write out the equation ( |96| ) and compare like derivatives. 
The collection of coefficients of the derivatives of the g k , starting from the highest deriva- 
tives, form a set of recursive relations for the coefficients in the intertwining transfor- 
mation. For example, one immediately sees from the d^gi and d 2 g<i coefficients that 
b<i{f) = and 0.2(7") = 0. Thus #2 does not appear in the physical even-parity gauge- 
invariant variable. As each successive intertwining coefficient is determined, the remain- 
ing relations are simplified. 

In anticipation of future developments, it proves convenient to rescale the remaining 
coefficients as follows 

b 4 (r) = iV 3 6 4 (r), (97) 
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ai(r) = A^ 3 ai(r), 
a 3 (r) = Na 3 (r), 
a 4 (r) = iVa 4 (r). 

It is also necessary to add zero in the form of an arbitrary multiple Ci(r), C2(r) of the 
(unsealed) momentum constraint variables c r , c g e minus their expressions in terms of 
the gk- This reflects the fact that the right-hand side of (^) may contain an arbitrary 
amount of momentum constraint variables. The Hamiltonian constraint variable is not 
added as it would mix in metric perturbations which are not present. [Note however 
that the Hamiltonian constraint variable appears on both sides of (|96|) because of its 
presence in ([75])- (|75D .] 

After these simplifications, reduce second and higher derivatives of N to first and 
no derivatives by taking derivatives of N 2 = 1 — 2M/r, and denote differentiation of 
coefficients with respect to r by primes. Then the set of recursive equations obtained 
from (|96|) becomes 

64 - rb' 4 (98) 
ra\ + ra 3 + b 4 

Xi(rai + ra 3 + 64) + 2ra 4 — 2r 2 a' 4 

+Arb 4 NN' - 8r%NN' - r 2 b'^N 2 - A^rci/2 
ra\ + 64 + r 2 a' 3 — rc\/2 

raiN 2 + 26 4 iV 2 - r 2 a[N 2 - 2r 2 fii^' - 2r~b A NN' + Nr 3 c 2 /2 

-ra 4 + 3ra 4 N 2 + 2ra 3 N 2 + 2b 4 N 2 

-2r 2 ~ ai NN' - 5r 2 a 3 NN' - 2rb 4 NN' - r 3 c 2 N - r Cl N 2 

r\ £ {h 4 + axN 2 + 2S 3 iV 2 - ra 3 NN' - Cl N 2 /2) 

-2A^ 2 (r5! - ra 4 + 26 4 iV 2 + 2rb 4 NN' - 2r%N' 2 ) 
-Ar 3 a' 4 NN' - r 3 N 2 a' 4 ' + r 3 c 2 N 

-2ra 4 - 2rS 3 iV 2 + 4b 4 N 2 + 2r 2 a 3 NN' - 4r~b 4 NN' 



o r gi : 

r 

-^-0 r g 4 : 



2N\ 
2N 5 



-d r gi : 



N 3 



-92 ■ 



N 3 



-94 ■ 



N 3 
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-Ar%NN' - r 3 a'^N 2 + r Cl N 2 - r 3 c 2 N 
N 3 ~ 

— gi : AN 2 (rh A - r^N 2 - 2b A N 2 ) + N^N^lOr 2 ^ + 28r6 4 - Sr 3 a[) 

-%r 2 N 2 N' 2 (r~ ai + 6 4 ) - 4r 2 AWa 4 - r 3 ^ 4 
+A £ r Cl iV 2 /2 + 3r 3 c 2 iV 3 + r 4 c 2 N 2 N' 

From the <9 2 g 4 coefficient, one has the relation 

b 4 = rd r b 4 (99) 

from which one concludes 6 4 = r. (An overall multiplicative constant can be taken to 
be unity without loss of generality since it represents a constant scale of the variable 
A g .) To cancel the d r g 2 term, one must set 

S 3 (r) = -Si(r) - I. (100) 

Cancelling the d r g^ and d r g\ terms determines c\ and c 2 to be 

d(r) = 2(1 + fii - r9 r ai) (101) 
c 2 (r) = ^- - 2rN' + a^N - 2rN') - rNa'^j . 

The coefficient of g 2 determines a 4 in terms of ai 

a 4 = 2N 2 -rA^AT' + a^iV 2 -rNN'). (102) 
After imposing this, the coefficient of <9 r g 4 gives the equation 
-X e + AN 2 -lQrNN' + a l (-X e + 2N 2 -8rNN') + ra' l (X e -2N 2 + 2rNN') = 0. (103) 
Inspection suggests the substitution 

ai = ar — 2 

which leads to 

-QMa + X e = 0. 

Thus, one has 

5i = ?57 - 2 - ( 104 ) 
6M v ; 
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After applying this relation, all the intertwining coefficients are fixed. The remaining 
equations implied by the gi, g% and g± coefficients are satisfied identically if the lapse 
takes the Schwarzschild form N 2 = 1 — 2M/r. These final relations are the consistency 
conditions of the over-determined system of equations. It is easy to see that if the lapse 
weren't special, e.g. Schwarzschild, the intertwining would not have been consistent. 
Assembling the above results gives the physical even-parity gauge-invariant variable 

A g = N(rN 2 d r g 4 + N 2 {^ - 2) 9l + (1 - ^)g 3 + (- - ^ + ^)gX (105) 

This variable satisfies the Regge- Wheeler equation 

\ N 2 

-d 2 + N 2 d r N 2 d r - V RW {r)\A g + — c r (106) 
N 2 ,2\ t AM N 2 2M A £ , 

+ — (-Z )Cg8 + —{ -T-JCfc = 0. 

r z 6 r r l r 6 

As with the odd-parity equation, there are constraint violating terms present to show 
how the theory evolves off-shell (perturbatively). On-shell this is the Regge- Wheeler 
equation. The gauge-invariant variable A g has even-parity and is constructed from 
time-derivatives of the extrinsic curvature. 
One can reach the Zerilli equation 



-d 2 + N 2 d r N 2 d r - V z (r))A g + • • • = 0, (107) 

where A g = DiA g and the Zerilli potential is 

2N 2 ( \ 

V z (r) = —, — - (n + l)n 2 r 3 + 3Mn 2 r 2 + 9M 2 nr + 9M 3 , 

r d (nr + 3M) 2 \ J 



108) 



(n = (£ — l)(£+2)/2) by applying the appropriate intertwining operator D\ ( |I4D to ( |106| ) 
and using (p^|) . Note the Zerilli potential has an unphysical singularity at nr + 3M = 0. 
Making this final transformation is unnecessary and arguably undesirable because of the 
unphysical singularity structure of the potential. Both the odd- and even-parity gauge- 
invariant variables satisfy a Regge- Wheeler equation, so parity is not the explanation for 
the existence of the Zerilli equation. Whether the existence of the Zerilli equation has a 
deeper meaning than accidental coincidence is an open question. From the standpoint 
of computation, however, the physics of even-parity perturbations is fully captured in 
the Regge- Wheeler equation without spurious unphysical complications, and the choice 
of description should be clear. 
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9 Conclusion 



For completeness, the full set of gauge- invariant perturbation equations (^), (|68|) , (jS5f), 
([86l), (H), ( P-06|) is gathered here 

N 2 \ 
-d 2 + N 2 d r N 2 d r - — {X e - QrNN')\A u - 2N 2 c u9 = 0. (109) 



N 2 \ 

-8 2 + N 2 d r N 2 8 r - — (A, - QrNN'))A g (110) 

, N 2 N 2 .2X t AM N 2 .2M X e 
+— c r + — (— )c 9 e + — ( — )c h = 0. 

ry £* ry ■£* -t ry ry£i ry ]q 



N 2 

-d 2 + N 2 d r N 2 d r - — X e )c ue = 0. (Ill) 



-d 2 + N*8 2 + (^ + N*N') dr ~^)c h+ *™c = 0. (112) 
iV 2 _ x 2 iV 4 A, X f N 3 N' s 



+ N 2 d r N 2 d r - ^ (2 + A, + 2rJVJV')) c r + ^^c 9 , + (WiVU - ^^-) c, = 

(113) 

\ J\j2 v OAT 

+ N 2 d r N 2 d r - - A ^)c g0 + -^{N + rN')c r - N 3 N'd r c h = 0. (114) 



ry 2 / iy^ 



Several important lessons are learned from the example of perturbing the Schwarzschild 
black hole. The one specific to the Schwarzschild case is that the Regge- Wheeler equa- 
tion suffices to describe the evolution of both odd and even parity physical perturbations. 
This result has been known mathematically for some time|l], [14], 122] , but a skeptic might 



assign significance to the observation that the Zerilli equation always seems to arise 
from even-parity perturbations [||, ||, [|. The calculation presented here should answer 
that doubt. The Regge- Wheeler equation has been obtained directly in even parity by 
an intertwining transformation, which decouples the equations and involves the fewest 
number of derivatives possible. The intertwining procedure is an effective tool which 
systematizes a technically involved computation and makes it transparent. The Zerilli 
equation can be reached by a further transformation which introduces an unphysical 
angular momentum-dependent singularity into the potential. Chandrasekhar[|I|] (pp. 
198-199) discusses whether one should dismiss the Zerilli equation and concludes not. 
"Dismiss" is probably too strong a word, but one certainly need never use the Zerilli 
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equation in computations, as all of the physics in it is captured more succinctly in the 
Regge- Wheeler equation. 

Perhaps the principal lesson which follows from the Schwarzschild example is the 
recognition that working with the fourth-order Einstein- Ricci formulation shifts the em- 
phasis from metric perturbations to curvature perturbations and in so doing reveals that 
the physical gauge-invariant quantities propagated by the Regge- Wheeler equation are 
formed from curvature perturbations and their spatial derivatives. Furthermore, one dis- 
covers that the gauge-invariant perturbations in unphysical directions are given by (time 
derivatives of) violations of the linearized Hamiltonian and momentum constraints. Ad- 
ditionally, the constraint violations evolve as a closed hyperbolic subsystem ( lll )- flni ), 
so that if the constraints are satisfied initially, they continue to be satisfied. It is worth 
emphasizing that this evolution subsystem for the constraints comes directly from the 
dynamical equations and not from separately considering the contracted Bianchi iden- 
tities. The information for evolving the constraints is embedded within the dynamical 
part of the theory, as it must be for consistency. 

These observations are general and help one to organize the calculation of gauge- 
invariant perturbations in a wide class of constrained hyperbolic theories: After recast- 
ing a constrained theory in hyperbolic form without gauge-fixing, say by the procedure 
outlined in [Q, one perturbs a complete set of independent combinations of the fun- 
damental variables which vanish in the background. Since the combinations vanish in 
the background, their perturbations are necessarily gauge- invariant |T5| . The constraints 
themselves provide a subset of these variables, and their perturbations will evolve as a 
closed subsystem if the equations are consistent. From a complete set of gauge- invariant 
variables, the physical subset which are independent of the constraint variables can be 
constructed in principle by using the intertwining procedure described in the text. 

The fact that the gauge- invariant perturbations of Minkowski space and the Schwarzschild 
black hole are curvature perturbations which split into the true linearized degrees of 
freedom and constraints encourages speculation about the nature of the true degrees 
of freedom in the fully nonlinear theory. From the twenty degrees of freedom of the 
Riemann tensor in four dimensions, the ten components of the Ricci tensor are obtained 
by tracing with the spatial 3-metric must of course be added to R k ikj to obtain 

the spatial Ricci tensor). Respectively, this accounts for one, three and all six degrees 
of freedom of Roioj, Roijk and Rijki- The ten Ricci components vanish (in vacuum) by 
equivalence of the Einstein-Ricci formulation with Einstein's theory. Six more degrees 
of freedom are obtained, three each, from divergences of Roioj and Ro^k- By the Bianchi 
identities, these reduce to combinations of derivatives of the Ricci tensor and hence also 
vanish. This leaves two degrees of freedom each in Roioj and Roijk, namely the transverse 
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traceless parts. Work is in progress to prove the conjecture that something related to 
the transverse traceless parts of the Riemann tensor are the true degrees of freedom of 
the gravitational field [P^]. 
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